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We update the single-coupling bounds on R-parity violating supersymmetry using the most up to 
date data as of October 2009. In addition to the data listed in the 2009 Review of Particle Properties 
[T], we utilize a new determination of the weak charge of cesium- 133 _2;, and preliminary r-decay 
branching fractions from Babar [31. Analysis of semileptonic D-decay is improved by the inclusion 
of experimentally measured form-factors into the calculation of the Standard Model predictions. 
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I. INTRODUCTION 



R-parity violating supersymmetry (SUSY) interactions |H IS] provide a convenient framework for quantifying quark- 
and lepton-flavor violating effects that new physics beyond the Standard Model (SM) may have, independently of 
whether SUSY truly exists in nature or not. Consequently, various authors have used a variety of flavor sensitive 
observables to constrain the sizes of these couplings [5J |SJ |71 |S1 IHl UHl HU ttH ESI IHl [151 US]. In this paper, we 
update the single-coupling bounds, namely the bounds on the individual R-parity violating couplings when only that 
particular couphng is assumed to be non-zero, using the most up to date precision data available as of October 2009. 
These include various lepton and meson decay ratios, CKM matrix elements, and the weak charge of atomic nuclei. 

The superpotential of R-partiy violating SUSY interactions is given by [U [5] 



(1) 



Here z, j, k are generation indices, while S'i7(2)-weak isospin and S'[/(3)-color indices are suppressed. The coefficients 
are antisymmetric in the first two indices, while A"^^. are antisymmetric in the latter two. Consequently, there 
are 9 independent LLE couplings, 27 independent LQD couplings, and 9 independent UDD couplings. The A^^-^ 
couplings lead to baryon number violating effects and are already very strongly constrained by proton decay [S], 
either individually |17] or in products with the Xijk and AJ^j, couplings [18 , so they will not be considered here. 
The explicit forms of the LLE and LQD interactions in terms of four-component spinors are 



C 



LLE 



= A 



C 



LQD 



(i<j)k 



(j>iLekB.ejL + e^LSkRViL + e.^B.^'iL^oi^ ~ (j^jLekLtfiiL + eiLekRMjL + Kb^^'^l^h^ 

{viLdkRdjL + d]LdkRViL + dln^fLdjL^ - [ciLdkRUjL + U^LdkR^iL + d*kR^iL^]l 



- h.c. 
+ h.c. (2) 



Note that the charge-conjugated chiral fermion fields are denoted — {Il/rY — {f'')R/L- The exchange of squarks 
or sleptons mediate interactions among the SM fermions, and the strength of these interactions will be proportional 
to the ratio of coupling constant squared to the exchanged sparticle mass squared. To simplify or notation, we follow 
Ref. i5 and define 



^ijk 



1 



|Aijfe|' 



'i]k 



1 



IKjkl'^ 

M| 



(3) 



Shifts in various observables will be expresses in terms of these dimensionless parameter combinations, 
bounds on the coupling constant will be shown with all the sparticle masses set to 100 GeV. 
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In the following sections, we look at the bounds from /i and r leptonic decays, t ^ nv and tt decays, CKM 
unitarity, semi-leptonic D and leptonic decays, and the weak charge of cesium-133. For analyses involving r- 
decay, the impact of preliminary r-decay data from Babar [3] is discussed. The analysis of semi-leptonic D decay is 
improved by a new calculation of the Standard Model (SM) predictions which include the effects of experimentally 
determined form-factors. The analysis of the weak charge of cesium-133 corrects an error in Ref. Bounds from 
Z-peak observables are not updated since no new data have been generated since the 2005 review by Barbier et al. 
[S]. The bound from neutrinoless double beta decay [B] will be discussed in a separate paper |41j . 



II. fi AND r DECAY 




FIG. 1: Possible R-parity violating contributions to (a) ^ 



'Ueiy^, (b) 



e VeVT, and (c) r — > VfiVr- 



The LLE couplings Xij^. affect the decays /i — s- e t'e^'^, r — > e i^e^^n E^nd r ^ /i f^^i^r- via the processes shown 
in Fig. [ij The operator induced by the exchange of CkR {k = 1,2,3) via the coupling A(i<j)fc is 



lA 



ii<j)k\ 
M? 



Wt<3)k\ 

2M? 



This will interfere with the SM operator 

_ 4G 
V2 

shifting the effective coupling to 



4G, 



4G, 



\/2 \/2 

In particular, \i2k will shift the muon decay constant to 



1 + ''(j<j)fe(efefl) 



Gn 



Gi 



1 + ri2k{ekR) 



and this shift will also affect other observables to be discussed later. The ratios 



and Rt 



will be shifted to 



1 + 2 



Rt — [Rt] 



SM 



+ 2|ri3fc(efcj^) - r23fe(efeK)| 



The SM predictions for the decay widths including radiative corrections [13 HO] are: 



SM 



m 



MM^ 967r3 



m 
m 



(4) 



(5) 



(6) 



(7) 



(8) 



(9) 
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phase space 


W propagator 


photon 


r(r- 


flU^ Vr) 


f (ml/ml) = 0.9726 


5w = 1-0003 


5; = 0.9957 


r(r- 




f (ml/ml) = 1.0000 










f {ml/ml) = 0.9998 


5'^- = 1.0000 


St; = 0.9958 



TABLE I: The corrections to the leptonic decay widths of the r and ^. 



V{t ev^Vr (7))sM 
r(/i ^ ev^v^ (7))sM 
in which /(x) is the phase space factor 



/ TO 



967r3 ' \ ^ ^ ^ 



2\ 1 



/(x) = 1 - 8a; + 8x3 - x" - 12x^1^^ ^ 



5^ is the W propagator correction 



w 



^ ^ 3 TO? 



(5^ is the radiative correction from photons 

^ 27r V 4 

and the values of the running QED couphng constant at the relevant energies are |20] 

a-i(TO^) = _ A In !!V + 1 ^ I36.0 , 
ATT rrif, DTT 

q;~^(to7-) « 133.3 . 

The numerical values of these corrections are shown in Table |lj The SM predictions for the ratios are therefore 



(10) 



(11) 



(12) 



(13) 



(14) 



[Rr ]SM = r(^^A^^A.'^r(7))sM ^ fi^l/^l) ^ ^ ^ ^ ggg ^ 

TM SM r(^->ePei/^(7))sM m5 f {mil ml) (5(^ 



SM 



r(T ^ ez/gi^r (7))sM _ fiml/ml) 



= 1.028 



r(r ^iv^Vr (7))sM f {ml/ml) 
The experimental values of these ratios from the Review of Particle Properties fT] are 



(15) 



■ B{t llV^Vr (7)) 



(2.197034 ± 0.000021) x lO^S 



Tr (290.6 ± 1.0) X 10-15 s 
6(t ^ ePei/^ (7)) _ (17.85 ±0.05)% 

B{t ^iv^vr (7)) ~ (17.36 ±0.05)% 



(17.36 ±0.05)% = (1.312 ± 0.006) X 10\ 
1.028 ±0.004. (16) 



The effect of a —13% correlation between B{t — s- ei/g 1/^(7)) and B{t — > ^v^Vt (^y) on the error on i?^ is small. 
Allowing only one of the A's to be non-zero at a time, comparison of Eqs. ([9|, (15 1, and (16 1 places the following 2(T 
bounds: 

lAi^^lf"^) < 0.05 [i?.,], 



|Al3fc| 



< 0.05 [i?.] 



|A23fc| ( ^^S^) < 0.05 [Rr] , 0.06 [Rr^] 



A new but still preliminary value of Rr from Babar was announced at ICHEP 2008 [3] as 

1 



[Rr_ 



Babar2008 



0.9796 ± 0.0038 



= 1.021 ±0.004 . 



(17) 



(18) 
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Including this value will change the world average to 

Rr = 1.025 ±0.003 , 

and the corresponding 2a bounds will be 

We see that the bound on Aia^ will be tightened. 



(19) 



III. TT AND r DECAY 



Ul 



— t— 



(a) 




Ul 




(b) 



FIG. 2: Possible R-parity violating contributions to 7r~ {£ = ei = e, or 62 — /i) that interfere with the SM amplitude. 

The indices are i = 1 or 2, fc = 1, 2, or 3 in (a); while (jk) = (11) or (22), with i = 3 — j or 3 in (b) due to the anti-symmetry of 
Xijk in the first two indices. The interference of (b) with the SM amplitude is suppressed due to the smallness of the electron 




(a) 




FIG. 3: Possible R-parity violating contributions to the decay t ^ tt fr. 

Possible R-parity violating contributions to the decay 7r~ {£ = e or fi) are shown in Fig. |2] and those to 

— > Tr~h'r are shown in Fig. [sj Since we are only interested in placing bounds on the individual R-parity violating 
couplings separately, we will ignore the (b) diagrams in both cases. 

The processes of Fig. [2]^ a) and Fig.jsj^a) induce the following operators: 



[!«): 



dkR 



Ficrz 



|Ailfe| 

'2MT 



^31k\ 



2Mi 



-{''fLl'"<L){uL7,.di 

R 



2M? 

dkR 
|-^31fcP / 



These interfere with the SM operators given by 



V2 



and 



V2 



(20) 
(21) 



and shift the 7r-decay widths to 



r(7r- iTvi:) = [r(^- ^ (.-w-)]^y, 

while the r-decay width is shifted by 



\Vud 



kR) 



r(r TT l^r) ^ [T{t -^n i^r)]gM 1 + |^ ^31fc ('^fei? ) 

Here, we have neglected any relative phase between the SM and RPV contributions. The ratios 

Rtt = 
Rttt — 



r(7r- - 


e 




r(7r- - 






r(r-- 






r(7r- - 







are shifted to 



Rtt — [Rt:]sM 
Rttt = [Rt7t]sM 



+ 1^ [r'nkidkB.) - 7-2uKfl)} 
1 + 1^ { i'2ik{dkR) ~ r'^ik{dkR) } 



At tree level, the SM prediction for the 7r-decay widths is given by 

2 2 

[r(7r- ^ri7«)]sM,trce - (V2G F\Vud\) 

while that of r-decay into tt Vt- is 

[r(T" ^ TT" I/^)]gM,troo = [^20 f\Vu. 

where the pion decay constant /^r is normalized as 

(0| u7m 75^(0) |7r"(q)) = ig^/, 

Taking ratios, we find 

[r(^- 



32^ 



1 - 



J TT ' 



J TT 



[Rtt ]sM,trcc 
[Rttt ]sM,tree 



e t^e)]sM,trcc _ ™| (1 - ml/ml)^ 

[T{tt- U^) ]sM,trcc ml (1 - ml/mlY 

[r(T^ ^ TT^ Z/,-) ]sM,troc 



(1 — ml/ml)^ 



[T{tt t';,)]sM,troc 2m2m^ (l-m^/m^)- 



1.283 X 10"^ , 
= 9.756 X 10^ 



Radiative corrections to these relations have been calculated in Ref. |^ and modify them to 

[-RttJsM — [^TT ]sM,trcc (1 + (^^tt) , 
[^tttIsM = [Rttt]sMMcc{^ + SRttt) , 



with 



SR^ = -0.0374 ± 0.0001 , (5i?^^ 
The uncertainty in these corrections is due to the uncertainty from strong interaction effects. Therefore, 



1.235 X 10" 

+0.009 iri3 



[Rtt]sM 

[Rttt]sm = 9.771+^;^^^ X 10 



On the other hand, the current experimental values are [T] 
S(7r^ePe(7)) (0.01230 ± 0.00004)% 



Rti 



B{tt fiD^, (7)) (99.98770 ± 0.00004)% 



(1.230 ± 0.004) X 10" 
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R-r 



B{t ^TTPr (7)) (2.6033 ± 0.0005) x lO'^s (10.91 ± 0.07)% 



Tr B{7T fi (7)) (290.6 ± 1.0) X 10-i5s (99.98770 ± 0.00004)% 



(9.775 ±0.071) X 10^ 



The magnitude of the CKM matrix element Vud is [22] , 

= 0.97418 ± 0.00027 . 



(33) 



(34) 



Comparison of Eqs. ( 25 1 , ( 32 1 and ( 33 ) leads to the following 2a bounds assuming only one of the couplings is non-zero 
at a time: 

lAUl f < 0.06 . 



Another preliminary result announced at ICHEP 2008 from Babar [3] was 

= 0.5945 ± 0.0063 



Babar2008 

Using the current world average value of B{t^ e^v^v.,-) — (17.85 ± 0.05)% [T] we find 

[S(r-^.-O]Baba.2008 = 10.61±0.12. 

Including this value will shift the world average to 

B{T-^T:-v-r) = 10.83 ±0.06, 

and the ratio i?^^ to 

R^^ = (9.703 ± 0.063) X 10^ . 



(35) 

(36) 

(37) 
(38) 
(39) 



The error will be reduced somewhat and the central value shifted down by about Icr. The 2a bounds will become 

(40) 



|Aiu.| < 0.04 [i?..] , |A^,,| f ^ ) < 0.08 [i?..] , 



My 



the change mostly due to the shift in the central value of Rrv 



IV. CKM UNITARITY 

The SM values of the CKM matrix elements Vud, Vus, and Vub must satisfy the imitarity relation 



w. 



SM,2_^|^SM|2 



W: 



SM|2 



= 1 



(41) 



Deviation of the measured values from this relation could be a sign of new physics. The value of \Vud\, cited above 
in Eq. ( 34 1 , is obtained from the comparison of superallowed 0+ 0+ nuclear beta decays and muon decay [22] , the 



former used to extract the product GpVud and the latter used to cancel the Gp- The SM operators relevant for these 
decays are 



4Gf 

71" 



and 



V2 



(42) 



Any new physics amplitude which interferes with these operators will affect the extracted value of V^d- The first 
operator is the same as the operator responsible for the decay 7r~ — > e~vZ. The RPV amplitudes which interfere with 
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FIG. 4: Possible R-parity violating contributions to nuclear beta decay. 




FIG. 5: Possible R-parity violating contributions to semileptonic if-decay. 



this are the same as those shown in Fig. [2] except with the ui^ hnes pointing toward the future as shown in Fig. |4] 
Since Fig. Qb) depends on two different RPV couphngs (its interference is also suppressed by the electron mass) we 
will only consider Fig. |4][a) for which the corresponding operator is 



M? v-L-iyv-eL-^y ■ ' 2M? 



2Af? 

dkR 



K7^rfL)(eZ7Mi^ei) . (43) 



We have already discussed muon decay in section |TT] where we found that Gp will be shifted by ri2k{ekR), cf. Eq. Q. 
Therefore, the shift in Vud will be 



SM|2 



ud 



V, 



ud\ 



AikidkR) - 2ri2fc(efcH) 



The values of V^s and V"„f, extracted from semi-leptonic K and B decays are 

\Vus\ = 0.2255 ±0.0019 , 
iKbl = (3.95 ±0.35) X 10^3 . 



(44) 



(45) 



The RPV diagrams that contribute to semi-leptonic iiT-decay are shown in Fig. [5] Similar diagrams contribute to 
semi-leptonic B-decay. None of these diagrams depend on a single RPV coupling so we may neglect them and assume 



V = 

US us 1 



V u ^ v 



SM 



(46) 



Then, 



1 



\Vud 



T'iik{dkR) - 2ri2fe(efcfl,) 



l-\Vus?-\Vub? 



The 2(7 bounds on the couplings are 

'lOOGeV' 



Ms 



< 0.03 [Vud] 



|Ai 



2fe 



/lOOGeV 



0.9999 ±0.0011 



< 0.03 [Vud] 



(47) 



(48) 



8 



V. SEMI-LEPTONIC D AND LEPTONIC D^-DECAY 



cl 




(a) 



Cl 



sl 



-t 



A/* 
i2k 

(b) 



ViL 



FIG. 6: Possible R-parity violating contributions to semileptonic D-decay and leptonic Ds-decay. 

The process shown in Fig. [6][a) affects semileptonic D-decay, while that in Fig.|6jb) affects leptonic il's-decay. They 
are both described by the same operator given by 



This interferes with the SM operator 



2Mi 



-{e'iLl''vtL){sLlt.CL) 



(49) 



(50) 



shifting the D and Ds decay widths by 



SM 



[T{D ^ K*£iUeJ] 



SM 



= l + , (51) 



SM 



Vr 



where we have neglected any relative phase between the SM and RPV contributions. Following Ref . jS] , we define the 
ratios 



Rd+ 









B{D^ - 


e+ 




B{D+ - 


-/^^ 




B{D+ - 






B{D+ - 






B{D+ - 




-i.ei^(892)0) 



the shifts of which are 



R 



D+ 



Rd+ 



[Rd' 



o\SM 



[R 



D + JSM 



\SM 



(52) 



(53) 



The experimental values of these ratios are currently [I] 

(3.32 ±0.13) X 10^2 
(3.61 ±0.05) X 10-2 

(9.4 ±0.8) X 10-2 
(8.50 ±0.26) X 10-2 
(5.4 ±0.5) X 10-2 



Rd" — 
Rd+ = 



R*D+ - 



= 0.92 ±0.04 
= 1.1 ±0.1 , 
= 0.98 ±0.10 



(54) 



(5.51 ±0.31) X 10-2 

Calculating the SM predictions of these ratios requires knowledge of the form-factors for the matrix elements [21] 

(K| 57*^(1 -75)c|i?) and (if* | S7^(l - 75)0 , (55) 
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for which good experimental data now exist from FOCUS |25l [27] , Belle 
Details of our calculation are presented in the appendix. The results are: 



], Babar US], and CLEO [301 III] ■ 



[Rdo]sm = [Rd+]sm = 1-04 ±0.02 (la), 1.04 

1 -1 



R 



D+ \SM 



+0.02 
'0.06 
-0.0 
-0.007 



1.060i° ™5 (2a) 



(56) 



The analysis of Ref. [S] used the value of (1.03) ^ without any errors for all three ratios. This would be the value 
of Ruo and Rd+ if form factors are ignored and the D and the K treated as point particles with the interaction 

{K^ dfj,D)W'^, and it lies within our calculated range above. For the ratio R*]j+ , if form factors are ignored and the D 
and the K* treated as point particles with the interaction K'^^DW^, its value would be (1.12)^^, which illustrates 
the importance of taking form-factors into account. The value of \Vcs\, also extracted from D semileptonic decays 
and Dg leptonic decays [25], is 



IK, I = 1.04 ±0.06 



(57) 



and for our current purpose we can set it to one. In comparing Eqs. (^3|, (54 1, and (56), we allow the SM predictions 



to scan the entire 20" range of Eq. ( 56 ) and pick up the weakest bounds on the couplings. The resulting 2a bounds 
are 

[^^^^) < i^Do] , 0.2 [Ro.] , 0.2 [R*^,] , 



22k\ 



100 GeV 



< 0.1 [Rdo] , 0.4 , 0.3 [R*a^] . 



(58) 



Next, define 



The shift of this ratio is 



Rd, (r/i) 



Rp, {ry) 

[i?C,(T^)]g^ 

The current experimental value is 



1 + 



B{D+ ^ /i+i.^) 



Vr, 



RoATfi) = 

while the tree-level SM prediction is 

m 

RdA'tN 



(6.6 ±0.5) X 10- 



(6.3 ±0.5) X 10-3 
2 (1 - mllm\jy 



Comparison of the two leads to the 2a bounds given by 

' 100 GeV ' 



= 10.5 ±1.1 



9.76 ±0.03 . 



lA' 



22k\ 



Ah 



100 GeV 



< 0.3 [i?i5,(TM)] 



(59) 



(60) 



(61) 



(62) 



(63) 



VI. ATOMIC PARITY VIOLATION 

The diagrams shown in Fig. ([t]) lead to the following effective couplings between the quarks and the electron: 

dkR dkn dkn 
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-T 

I 
I 

I 
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(a) 



ul 





(b) 



FIG. 7: Possible R-parity violating contributions to atomic parity violation. 



(&) 



eLdR) {dneL) 



2M? 



The parity violating parts of these interactions are 

IA' P r 

(a) : [(67^756) (u-f'^u) + (e7^e) {u-f^-fsu 



The parity violating SM interactions, on the other hand, are 

[C'i9(e7M75e)(g7^(?) + C2g(e7^e)(g7^75g) 

q—u,d 

with 



Clu — 
Cld = 



1 4 2 

2 + 3-^ ' 
1 2 , 



C2u = -^ + 2s' 
Cad - +^-2s' 



(64) 



(65) 



(66) 



(67) 



at tree-level, being the shorthand for sin^ 6w The weak charge of an atomic nucleus with Z protons and N 
neutrons is given by 



{Z.N) 



Ciu{'2.Z + N) + Cid{Z + 2N) = (l-4s2)Z-7V. 



(68) 



Note that (2Z + N) and {Z + 2A^) are respectively the number of up and down quarks in the nucleus. Since the 
presence of the above R-parity violating couplings will shift the Ci couplings to 



Ciu Ciu — r'l-j^kidkR) , 
Old Cid + r[ji{ujL) , 

the weak charge will be shifted by 

SQw{Z,N) = +2\{2Z + N)r[,^{dkR)-{Z + 2N)r[^,{ujL) 



(69) 



(70) 



Furthermore, since the quantity that is actually measured in atomic parity violation (APV) experiments is the product 
GpQwi Sind Qw is extracted by setting Gp equal to the muon decay constant G^, the LLE coupling ri2kiekR.) can 
also affect Qw via Eq. ^ : 



GfQw r\ \^ (~ \^ 

— ^ = Qwy^ - ri2k[ekR) \ 



(71) 
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For cesium- 133, the total shift will be 

(5Qh^(55, 78) = 376 r[^kidkR) - 422 r[jAujL) - Qw{55, 78) ri2fe(gfefl) . (72) 

(This formula differs from that provided on page 82 of Ref. [5 which contains a typo: the factor of —2 on the right- 
hand-side should not be there. This error seems to have propagated into the bounds listed in Eq. (6.47) on the same 
page. We correct for this error in quoting the bounds from Ref. [5 in Table [ll|) The latest experimental value of the 
weak charge of cesium- 133 is [2 

QwCUCs) ^ -73.16 ±0.29 ±0.20 = -73.16 ±0.35. (73) 

On the other hand, the SM value provided in the Review of Particle Properties [T] is 

[Qv^(iiCs)]3^ = -73.16 ±0.03, (74) 

which is based on a global fit to all electroweak observables with radiative corrections only from within the SM. 
Assuming further radiative corrections from new physics is negligibly small, and saturating the difference between the 
two with the RPV contributions, we find the following 2a bounds: 

, /lOOGeV\ , , , /100GeV\ , , /100GeV\ , , 

A'n, -jrr < 0.04 [Qw] , X[,i -rr^ < 0.03 [Qw] , Xi2k -jr < 0.08 [Qw] ■ (75) 

VII. SUMMARY 



Coupling 


Observable 


Previous 2a bound [5] 


New 2a bound 


with Babar preliminary T 






0.05 


0.03 








0.07 


0.05 






Qw{Cs) 


0.11* 


0.08 




Al3fe 


Rt 


0.07 


0.05 


0.03 


A23fe 


Rr 


0.07 


0.05 


0.05 




Rth 


0.07 


0.06 






QwiCs) 


0.03* 


0.03 




Allfc 




0.02 


0.03 






R-K 


0.03 


0.03 






QwiCs) 


0.05* 


0.04 




\' 

^12fe 


RdO 


0.271" 


0.2 






Rd+ 


0.44+ 


0.2 






R*D+ 


0.23+ 


0.2 




-^21*: 


R-K 


0.06 


0.06 






Rttv 


0.08 


0.07 


0.04 


\' 

^22fc 


Rd" 


0.21+ 


0.1 






Rd+ 


0.61+ 


0.4 






R*D+ 


0.38+ 


0.3 






-Rds(tm) 


0.65 


0.2 




Asife 


Rttv 


0.12 


0.06 


0.08 


^32k 


-Rds('^m) 


0.52 


0.3 





TABLE II: Current single-coupling bounds on various RPV couplings compared with those given in Table 6.1 of Ref. [5]. AH 
sparticle masses have been set to 100 GeV. The values with an asterisk have been corrected for an error in Ref. [5]. Those with 
a dagger were calculated with the SM prediction fixed to (1.03)^^ without any uncertainties. 

In this paper, we have looked at a variety of single-coupling bounds that can be imposed on R-parity violation from 
particle decay ratios and atomic parity violation. Our results are summarized in Table |TTj Compared to the bounds 
complied in 2004 by Barbier et al. [5], the bounds have steadily improved. For the bounds from observables involving 
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T-decay, the improvement is due to the results of ALEPH f32] released in 2005. The bounds from observables involving 
D and Dg decay have improved due to new data from CLEO [30l ESI EH ES] and Belle |28(|36i- The new bounds from 
the weak charge of cesium-133 is due to the reduction of theoretical uncertainty in its extraction from experimental 
data 12 . 

Further improvements on the bounds from observables involving r-decay is expected as analyses of Belle and Babar 
data get under way, each with hundreds of millions of r+r^ pairs, as is evidenced by the effect of the preliminary Babar 
data from 2008 0] . The bounds from observables involving 7r-decay can also be expected to improve considerably in 
the near future as the PIENU experiment at TRIUMF [3Z, and the PEN experiment and PSI [38 start their physics 
runs this year, eventually improving the error on i?^ by almost an order of magnitude. 

The bound on X'm from neutrinoless double beta decay will be updated in a subsequent paper |41j . 
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APPENDIX A: SEMILEPTONIC D DECAY RATIOS WITH FORM-FACTORS 



The semileptonic decays D Kiv {D+ K"e+iye or D° 
D° K*{892)-£+Vi) proceed via the SM interaction 



C 



V2 



K-t+vi) and D KHv {D+ K {m2fi+vi or 
KU^7m(1-75)c][^?7'^(1-75K] . (Al) 



For both cases, the hadronic matrix element is expressed in terms of multiple form factors. In this appendix, we 
calculate the ratios Rj^o , , and , defined in Eq. ( 52 ) , taking the latest experimental data on these form 
factors into account. 



D -> Klv 



For the decay D Kiiy, the form factors are defined as [53] 

{K\sj^{l - j,)c\D) = (Pd + Pk)^. f+{t) + {Pd ~ Pk)^ f~{t) , 
where t — {Pd — PkY ■ The decay width in terms of these form factors is 



(A2) 



247r3 



where we have defined 



2m|) 



1 



t 



2t 



\f+it)\' 



3 mj ( Ai/2(t,m|),m^) 



2ml, 



Mt) EE f+it) + —^ ^Mt), 

rajj — rrii 



t ( ml 



(A3) 



(A4) 



and 



A(a, b, c) = + 6^ + - 2ab - 2bc - 2ca . 



(A5) 



Since ^ to|,, the scalar form factor term is completely negligible for the electron, and the vector form factor 
can be extracted from the i-dependence of dT/dt. The result is fit by the single-pole function 



/+(o) 



(A6) 
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where /+(0) and mpoic are adjustable parameters. This form is motivated by the vector meson dominance model 
which prescribes the following expressions for the form factors |40] : 



m 



1 



cv 



Cv + cs 



1 - t/m% ' 



(A7) 



Here, cy and eg are constants, and mv and ms are the masses of the lowest lying D^-mesons with = 1 and 0+, 
respectively. Another popular form used to fit the /+ {t) data is the so-called modified-pole function 



U{t) 



MO) 



(1 - t/m%^){l -at/m\,,) 



(A8) 



where niD* = 2112 MeV, and /+(0) and a are adjustable. The fit values of mpoic and a from recent experiments are 
hsted in Table Iml 



Collaboration 


Decay Mode 


m-poio (GeV) 


a 


./-(0)//+(0) 


Reference 


FOCUS 




1.93 ±0.05 ±0.03 


0.28 ±0.08 ±0.07 


-1.7tii±0.3 


[25] 2005 


Belle 


D° ^K'i+ue {£ = e,ii) 


1.82 ±0.04 ±0.02 


0.52 ±0.08 ±0.06 




.28, 2006 


Babar 




1.884 ±0.012 ±0.016 


0.377 ±0.023 ±0.031 




.29, 2007 


CLEO 


D° K-e+Ve, D+ ^ K"e+iye 


1.93 ±0.02 ±0.01 


0.30 ±0.03 ±0.01 




Mi 2009 



TABLE III: Recent measurements of the D K form factors. 



The forms of the single-pole and modified-pole functions for the CLEO |31j central values of rripoic and a are shown 
in Fig. [s] They deviate from each other considerably beyond t ~ 2 GeV^, where data points are absent due to phase- 
space suppression. However, due to this same phase-space suppression, it turns out that the Rjj ratios are insensitive 
to this difference in the shape of f+{t) in this region of t, so we adopt the single-pole form for our purpose. 




FIG. 8: The single-pole function with mpoie 
end-point is at t = mjj ~ 3.5 GeV^. 



1.93 GeV (solid), and the modified-pole function with a = 0.30 (dashed). The 



FOCUS [25 , which looked at the muon channel, has determined the ratio /_(0)//+(0), in addition to mpoic and 
a, by performing a two-dimensional fit to cPT/dtdcos0(, where 9i is the angle between the neutrino and the kaon in 
the rest- frame. If we adopt the meson-dominance form given in Eq. (A7l for fo{t), with cy = /+(0), this ratio 
determines — cs/ cy as 



cs_ 

Cv 



polo 



-^l /+(o) 



(A9) 



In the FOCUS analysis, it was assumed that the ratio f-{t) / f+{t) was essentially independent of t, which would 
require ms ~ Wpoic. However, as noted above, data points do not extend into the regions of t which are phase- 
space suppressed, thus ms can be differ significantly from mpoio while still maintaining an almost constant ratio 
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f-{t)/f+{t) in the regions measured. Indeed, in Fig. ^we show how this ratio varies with t for the two choices 



TOpoio = 1930 MeV and and ms = m(i:)*o(2317) 



ms 

t < 2 GeV^, the ratio only varies from —1.7 to about 
to vary between mpoio and m(£)*Q(2317)*). 



2317MeV when /_(0)//+(0) = -1.7. In the region 



-1.5 even for the latter case. Thus, in our calculations, we allow 



0.0 



-0.5 



-1.0 



-1.5 




-2.01 



2 

t (GeV^) 



FIG. 9: The ratio f-{t)/f+{t) with ms = mpoic = 1930 MeV (solid), and ms = m(D*o(2317)=^) = 2317MeV (dashed). 
FOCUS central value is used for the value of ratio at t = 0. 



The 



The form-factors we use are therefore, 



Ml 

/-f(O) 

Mt) 
/+(o) 



1 



1 , v™l 



1 



</m| 



(AlO) 



For TTipoic we adopt the FOCUS value of 1.93 GeV, which coincides with the CLEO value. The errors on mpoic are 
small enough as to have no effect on Rd since they only change the form of f+{t) near the pole. The value of ro will 
depend on the value of the ratio /_(0)//4.(0) and our choice of ms as shown in Fig. 10 




/_(0)/A(0) 



FIG. 10: The dependence of ro on the ratio /_(0)//+(0) with ms = mpoic = 1930 MeV (solid), and ms = m{D*o{2317)'^) = 
2317 MeV (dashed). The shaded bands indicate the one- and two-sigma regions of /_(0)//+(0) as measured by FOCUS |25) . 



Substituting into Eq. ( A3 1 , we obtain the values shown in Fig. 
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for R 



D 



off from the graph, we assign the following la and 2a ranges to Rjjo and Ro+'- 

(i?Do)-i = (i?m)^' = 1.04 ±0.02 (la) , 1.04i°:°2 (2a) 



F(D Kev)/T{D K^iv). Reading 

(All) 
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1.08 



1.06 




/_(0)/A(0) 

FIG. 11: The dependence of (Rd)'^ = r(D Kev)/V{D Kfiv) on the ratio /_(0)//+(0) with ms = mpoic = 1930 MeV 
(solid), and ms = m(-Dso(2317) ) = 2317MeV (dashed). The shaded bands indicate the one- and two-sigma regions of 
/_(0)//+(0) as measured by FOCUS [25]. 



D KHv 



Next, we consider the decay D — > K*tvi. The form-factors are defined as 



{K*{PK,e)\-s^^{l-^^)c\D{PD)) 



c „ . PT c*^ 



mu + rriK 
+i |e* (niD + rriK') Ai{q^) - 



V{q^) 



{PD + PK)c.A2{q^)-^-^^qMmK')As{q^) 



q 



(A12) 



where 



2mK'A3{q^) = {ruD + mK*) Ai{q^) ~ {mo - rriK') A2iq^) 



AsiO) = Aq{0) 



(A13) 



and = (Pd — Pk)°' ■ Here, V{q^) is the contribution of vector intermediate states, while Ai{q^) and ^2(5^) 
are contributions of axial-vector intermediate states. A^^q^) combines with Ai{q^) and ^2(9^) to maintain the 
transversality of these axial-vector contributions. ^0(9^) is the contribution of scalar intermediate states. The 
condition ^3(0) = Ao(0) is necessary for the cancellation of the poles at q^ — 0. If we rewrite A^ in terms of Ai and 
A2, the above expression becomes 



(if*(PK,e)|s7a(l-75)c|i?(Pi3)) 



e- , p/3 p7 ,*S 



V{q^) 



+i |(m_D + riiK') 



(e* • q) qa 



Al{q') 



{^*-q) 

rriD + rriK' 



(Pd + Pk)c 



(to|, - m|-, ) qg 
g2 



A2{q') 



+i 



qai2mK*) Ao{q ) , 



(A14) 
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making the transversal! ty of the axial- vector terms manifest. The decay width in terms of these form-factors is 



r 



277r3 



dt 



3{mD + mK*Y 



^2 



A(m|,,m|-.,t) 



1 + 



3m|. 



2i 



T 



2i 

Al(t)A2(i) 



1 - 



+ 



6(m£) -t- TOif » )2m|- 
A3/2(m25,m2^.,t) 7 



t 



„2N 2 



t 



AM 



(A15) 



where t = q^. The Ao(t) contribution, being proportional to m|, cannot be measured at current experimental 
sensitivities. For the other form-factors, the following expressions motivated by vector-meson-dominance [39 is 
assumed: 



V{t) = ViO) 



1-t/ml ' 



(A16) 



The pole masses are taken to be my — m{D'*^) = 2.1 GeV and — to(Z?si(2460)*) — 2.5 GeV, which are the 
masses of the lowest lying Z?s-mesons with = 1^ and 1+, respectively. Experimental data is then fit with the 
ratios = V{0)/Ai{0), and r2 = A2{0) /Ai{0). These ratios have been measured by various experiments via the 
decay sequence D+ K* (892)° e+i^i, K*{892)° if-7r+, most recently by FOCUS |26], which also measured the 
decay sequence i4r*(892)^/i+f^, X*(892)^ K'^tt^ [27J. The current world averages are [T] 



r2 = 



V{0) 
Alio) 
A2{0) 



= 1.62 ±0.08 , 



= 0.83 ±0.05 



^i(O) 

For our purpose, we also need to specify Ao(t) which we assume is of the form 



Ao{t) = A3(0)+cs 



1 - t/ml 



with ms = rnj-i± = 2.0 GeV, the mass of the lowest lying D^-meson with 

Ao{t) _ , t/ml 



Ai{0) l-t/rng 



with 



(A17) 
(A18) 

0^. Normalizing to Ai(0), we have 

(A19) 



2mK* 



2mK' 



r2 = 1.09 ±0.03 



(A20) 



For ro, we arbitrarily assume that its Icr range is rg = ± rs. Substituting Eqs. (A16l and (A19) into Eq. (A15I, 
we obtain the values shown in Fig. 12 for the ratio — V{D — > K*ev)/V{D — > K* jiv). Reading off from the 

graph, we assign the following Icr and 2(T ranges: 



(A21) 
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1.050^ 
-3 



-2 



-1 





'"0 



FIG. 12: The dependence of = r{D K* ev) /V{D K* ^v) on the choice of ro. The wide vertical bands are our 

assumed one- and two-sigma regions of ro. The narrow curved bands indicate the one and two-sigma uncertainties due to the 
error in r2- The uncertainty due to errors in r^, mo-, and mx* are neghgible. 
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